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Abstract 

We study the recently proposed Covariant Renormalizable Gravity (CRG), which 
aims to provide a generally covariant ultraviolet completion of general relativity. We 
obtain a space-time decomposed form — an Arnowitt-Deser-Misner (ADM) repre- 
sentation — of the CRG action. The action is found to contain time derivatives 
of the gravitational fields up to fourth order. Some ways to reduce the order of 
these time derivatives are considered. The resulting action is analyzed using the 
Hamiltonian formalism, which was originally adapted for constrained theories by 
Dirac. It is shown that the theory has a consistent set of constraints. It is, however, 
found that the theory exhibits four propagating physical degrees of freedom. This 
is one degree of freedom more than in Hofava-Lifshitz (HL) gravity and two more 
propagating modes than in general relativity. One extra physical degree of free- 
dom has its origin in the higher order nature of the CRG action. The other extra 
qh| propagating mode is a consequence of a project ability condition similarly as in HL 

^ ' gravity. Some additional gauge symmetry may need to be introduced in order to 

get rid of the extra gravitational degrees of freedom. 

^ I PACS: 04.50.Kd (Modified theories of gravity), 04.60.-m (Quantum gravity), 

■ ll.lO.Ef (Lagrangian and Hamiltonian approach), 98.80.Cq (Particle-theory and 

. field-theory models of the early Universe) 

00 

1 Introduction 

In the recent years modified theories of gravity have attracted a considerable amount 
of attention. These modifications of General Relativity (GR) aim to improve the be- 
^ , haviour of the theory either at high energies (renormalizability, early-time universe) or 
at large distances (cosmology, galaxies), and sometimes at both regimes (see [H, 0] for 
such an attempt), in order to achieve a better agreement with observational data or a 
more plausible theoretical framework. Recently the so-called Covariant Renormalizable 
Gravity (CRG) was proposed in Refs. [3|-l5|. For a review of CRG and its comparison 
with other modified gravities, one can see p|. CRG aims to provide a power-counting 
renormalizable field theory of gravity that is covariant under spacetime diffeomorphism 
and possesses local Lorentz invariance at the fundamental level. Lorentz invariance of the 
graviton propagator of CRG is, however, broken dynamically at high energies. This is 
achieved by introducing an exotic fluid of unknown origin, which is coupled to spacetime 
in a rather complicated way. Horava-Lifshitz (HL) gravity 0] — another power-counting 
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renormalizable field theory of gravity — is based on tlie idea tliat space and time scale 
differently at high energies, 

x^bx, t-^¥t, (1.1) 

with a dynamic critical exponent z = 1, 2, 3, ... . This enables one to modify the ultra- 
violet behaviour of the graviton propagator to |A;|~^^, where k is the spatial momentum. 
In D spatial dimensions, z = D ensures the theory is power-counting renormalizable. 
Such spacetime admits a preferred foliation into spatial hypersurfaces, and hence the lo- 
cal Lorentz invariance is broken. Note that more general theories of such sort have been 
proposed in Ref. j^. CRG achieves a similar ultraviolet behaviour of the graviton prop- 
agator as HL gravity, but without introducing explicitly Lorentz noninvariant terms into 
the action. However, we should note that the renormalizability of CRG, as well as of the 
HL model, is assumed only based on the power-counting arguments. There are several 
potential pathologies that could ruin the renormalizability of this theory, such as gradient 
instabilities, ghosts or strong coupling. Since a violation of Lorentz invariance has never 
been observed, one could argue that CRG is a more natural modification of GR than the 
explicitly Lorentz noninvariant ones, in particular HL gravity and its generalizations (see 
e.g. [il,S,ai9|). 

Hamiltonian formalism provides a powerful tool for the analysis of constrained systems 
such as gravity. For example, it has been shown that HL gravity can be physically 
consistent {N > ) at high energies only if the projectability condition is imposed on the 
lapse function N [lO|: N = N{t). A similar result has also been obtained for the more 
general modified F(R) HL gravity [Uj. In this paper we seek to understand CRG from 
the point of view of Hamiltonian analysis. 

First we obtain the Arnowitt-Deser-Misner (ADM) representations of the two most 
interesting actions of CRG theory in Sec. [2] — one supposedly power-counting renormal- 
izable and the other power-counting super-renormalizable. Since the actions turn out to 
contain time derivatives up to fourth order, we consider ways to simplify the actions via 
introduction of additional scalar fields and partial gauge fixing. In Sec. [3] we analyze the 
actions using Hamiltonian formalism. Conclusions and some further discussions are given 
in Sec. H] 



2 ADM representation of covariant renormalizable 
gravity 

2.1 Action 

Let us review the action of the CRG theory proposed in jsj-Il] • In 4-dimensional spacetime 
the action of the supposedly power-counting renormalizable gravity is (corresponds to 
z = 3) 

Sg^s = j d^x^g 1^ - a {T^'^R.u + (3TR) (T'^^V^V. + iTVV p) {T^" R,u + /3Ti?)| . 

(2.1) 

We also consider the supposedly power-counting super-renormalizable gravity with the 
action (corresponds to z = 4) 

Sg,, = j ^'^v^ 1^ - « [(^^'^'^V^V. + iTWp) (T^^^R^u + /3Ti?)]'| . (2.2) 
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Here T^^, is the energy-momentum tensor of a perfect non-relativistic fluid with an equation 
of state parameter w 7^ 1/3, —1. We denote T = g^^^T^^. The actions for higher z could be 
analyzed similarly as these two interesting representatives of the odd and even z actions. 

In Ref. H it was proposed that the required perfect fluid can be realized with a 
constrained scalar field with the action 

S^ = j d'x^ I -A + f/(0)^ I , (2.3) 

where A is a Lagrange multiplier field. The energy-momentum tensor T^^ is defined to be 
that of an isolated scalar field in a potential V{(f)), 

Tt. = d,<pdA - g,u (^dp<pdPct> + V{cp)^ , (2.4) 

regardless of the form of the CRG action. The constraint implied by the action ( 12. 3p . 

^a^0a> + f/(0) = O, (2.5) 

will be imposed from the start. Thus the constraint (12. 5p . and consequently the action 
(12. 3p . vanish in the total action of CRG. According to the constraint (12. Sp . if we assume 
?7(0) > 0, the vector 5^0 is timelike. Then at least locally one can choose the direction 
of time to be parallel to (9^0, so that the constraint (12.51) yields 

Then on the flat background metric one obtains the energy density = U{(j)) + V{(j)), 
the pressure = f/(0) — ^(0) and the equation of state parameter w = p<f,/ associated 
with T^j^. For simplicity, it is assumed that V^(0) and ?7(0) are constants: 

f/(0) = f/o, V{<P) = Vo. (2.7) 

Now the action (12. ip can be written as 

5-3,3 = j d'x^ - " (9>9^0i?,,, + UoR) {d^<pd^<pV^V. + 2f/oV''Vp) 

X (a^0a>i?^, + [/oi?)} . (2.8) 

Likewise, the action for (12. 2 p is written as 

5-3,4 = j ^'^v^ - « W^d^'^V.V, + 2f/oV^Vp) (9>9>i?^, + UoR)f 

' (2.9) 

Here the constraint fl2.5p has been used and the following parameters have been deter- 
mined: 

^ 2(3w - 1) 2f/o - 4K) ' ^ 3u; - 1 2Uq-AVq' ^ ' ' 

In order to eliminate Uq from the actions (12. 8p and (12. 9 p we can use the constraint (12. 5p 
to write jsf 

d^4>d''(PR,, + UoR = d^<pd''<p (^R,, - ^g,uR^ , (2.11) 
9''0a>V^V, + 2f/oV^Vp = d^<j)d''<f) (V^V, - g^^^WVp) . (2.12) 
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Since (12. 3p vanishes, the total action {Sz = Sg^z + S^) for (12. 8 p can be written 



xd'^cpd'cp (v^V. - g^^uWp) d'^cpd^cp ( i?.A - ) • (2.13) 



For (12. 9p the total action can be written 



2n 



1 



2^ 



(2.14) 

The corresponding Lagrangians are denoted by L^, such that = J dtLz- 
2.2 Space-time decomposition 

We consider the ADM decomposition of the gravitational field (for reviews and math- 
ematical background, see [isl, |l3])- Assume that the spacetime admits a foliation into 
t = constant hypersurfaces E^, whose future-directed unit normal is denoted by n. The 
metric tensor ^%ab of spacetime is decomposed with respect to S^: 

gab = ^'^^gab - (^naUb , Tlan" = 6 , J = 1 , 

N = enat\ N'' = g\t\ = Nn'' + N\ (2.15) 

where e = — 1 (e = 1) for a spacelike (timelike) Sj, and the orthogonal projector from 
TM to TSt is 

g\ = S^b-^ri''rib. (2.16) 

The extension of tensors on to tensors on A4 is induced by the embedding of to A^, 
and by the projector map (I2.16p . Such extended quantities are, for example, the metric 
tensor gab of and the shift vector A^". From now on quantities defined on the spacetime 
Ai and associated with its metric ^^^gab are marked with the prefix We denote the 
covariant derivatives on Ai and by V and D, respectively. 

Taking the scalar field t as the time coordinate, when is the "flow of time" vector, 
and assuming is spacelike (e = —1), we can write the unit normal in terms of the lapse 
A^ and the shift vector A^*: 



rir 



/ 1 ]\fi\ 

-NV,t = i-N, 0, 0, 0) , n'^= (n°, n^) = (-, -—] . (2.17) 

The components of the metric of spacetime are decomposed 

^%oo = -N' + N,N^ , = ^%,o = N, , = g,j (2.18) 
and the components of its inverse are 

(4)^00 = -l/N^ , ^ iyo ^ ^^|^2 ^ (4)^ij _ _ (jy^VAT^) . (2.19) 

The natural volume element decomposes 

d^x^-^^)g = dtd^x^N . (2.20) 
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We assume c = 1 units. 

A 4- vector can be decomposed into components perpendicular and parallel to n'^: 

u>' = ^uf" + e(n ■ u)n>' , ^m'' = = - e(n ■ u)n^' . (2.21) 

For = d^(j) = by using fl2.17p we get 

= ^V'^^ = (0, D'cf) , Z^V = 9''D,<f) = g'^d.cj) , (2.22) 

and 

n-u = Vn(f) = n''V^^. (2.23) 
Note that for a vector jm'^ € TS^ and a covector jt"^ G T*Si we have the zero- components 

±M° = 0, ±t;o = iV*±^^^ • (2.24) 

This generalizes straightforwardly for tensors of any type defined on Ej. 

We will apply the ADM decomposition to the actions (12.131) and (I2.14p . The scalar 
curvature ^'^''R of spacetime decomposes as usual: 

= KijK'^ - + R + 2V^ [n^Wy - TfVy) 

= KijK'^ -K^ + R + 2V^ (n^'K) - ^D'DiN . (2.25) 

Here R is the intrinsic scalar curvature of and Kij is the extrinsic curvature of S^: 

= ^ {9^3 - '^D^iN,)) , K = g'^Kij , (2.26) 

where the dot denotes the derivative with respect to time t. 

The covariant derivative of the normal covector can be written 

V^n, = K^, - n^D, In N , (2.27) 

and it implies the following relations: 

V„n^ = D^lnAr, n'V^n, = 0, Vy = K. (2.28) 

These relations have been used in fl2.25l) . Note that we can also write 

V^{n^K) = K^ + VnK, (2.29) 

for example, in fl2.25p . 

The Ricci tensor ^"^^R^i^ of spacetime can be decomposed as 

^ ^Rfiu ~ ±Riiu ±Riin^i' ±Rni/^ii ~l" ^Rnn'^fi'^u i (2.30) 

where we have defined 



^^R^, = g'yj^^. = R^u + KK^, - 2K^pKP^ - ^^D^D^N + ^^CnuK^, , (2.31) 



^^iR^n = g'y^'^Rp, = D.K"^ - D,K , (2.32) 
^'2Rnu = 9'y^''%p = DpK"^ - D,K , (2.33) 
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- KijK'^ - (n^'K) + ^,DWiN . (2.34) 



In Eq. f l2.3ip R^u is the Ricci tensor of the hypersurface and Cnu denotes the Lie 
derivative along iVn^. Note that for any tensor field T that is tangent to Sj, £Ar„T is also 
tangent to S^. In Eqs. fl2.3ip - fl2.34p we have used the Gauss relation, the Ricci equation 
and the Codazzi relation, and in fl2.34p the decomposition of ^'^^R from fl2.25p was also 
used. Hence the Einstein tensor can be decomposed: 

= - = R^u + KK^, - 2K,pKP^ - j^D.D^N + ^CnuK,, 

- \g,u {r + K,,K''^ +K' + 2VnK - ^DW,N 
+ {D,K - DpKP^) + {D,K - D^K^^) + ]^n^n, [K^ - K,,K'^ + R) , (2.35) 
where we have also used f l2.29p . Thus in the actions fl2.13p and (12.141) we have 



R,, + KK,^ - 2K,,K''^ - ^D,D,N + ^CnuK.,, 



- ^g,, (^R + KkiK'^' + K^ + 2VnK - ^D'^D^N^ 
+ 2 {Vn<P) D'<P {D,K - D^K,,) + ^ (V„0)' {K^ - K,,K'^ + R) . (2.36) 

Let us then consider the covariant derivatives in the actions fl2.13p and (12.140 . The 
covariant derivative D of a {k, /)-tensor field T on is given by 

n rpui -Uk _ „0- I/l _ _ _ i^fe J3l _ _ _ A ^ rpai-ak /r, oyx 

-^M^ pi-pi - 9 ^,9 ai 9 ak9 pi 9 pi^aJ- ' 

where in the right-hand side one considers the extension of T on A^. For a scalar field /, 
Dfif = g^^uf implies the decomposition 

VJ = DJ-n^Vnf. (2.38) 

For the second order covariant derivative we obtain 

D,DJ = g'^g^VpV.f + g\^g%) [V pg\) Va/ • (2.39) 

Note that both D and V are torsion-free, which is the reason why the last term is sym- 
metrized over /i and v. Introducing fl2.16p into fl2.39p yields 

V^V,/ = D^DJ - 2n(^|V„V|,)/ - n^n.nPn'^V J 

- (V(^n^) + n,^p\Vnn\y) + n^^n''Vy)np + npUynPVnnp) Vnf ■ (2.40) 

Then we shall use the relations fl2.27p and (12.281) to obtain 

VpV,/ = D^DJ - KpyWnf - 2^(^(7'',) V„Vp/ + n^n,n''n-VpV./ . (2.41) 
As the last step in the decomposition we write 

V„ Vp/ = VpV„/ - V^VJ = VpV„/ - (ir/ - ripD"^ In N) VJ , (2.42) 
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which imphes 

g'J^nVpf = D^Vnf - K/Vpf , (2.43) 
n^n'^VpVJ = V„V„/ - {D' In N)DJ . (2.44) 

Thus the final decomposition is 

V^V,/ = D^DJ - K^^Wnf - 2n(^Z},)V„/ - 2n(^K^/V,/ 

+ n,n, (V„V„/ - [D' In N)D,f) . (2.45) 

The contraction of fl2.45p with (^^(yi'^^ gives 

V^/ = DW,f - KVnf - VnVnf + {D' In iV) A/ , (2.46) 
where we have also used 

^^^''D^DJ = g'W^D^f = D'DJ , ^^^''K^, = g'^K,, = K . (2.47) 
Then in the action we have 

- {Kij - gijK) Vnf + gij^n^nf] 

- 2 (Vn0) (A V. + K,,D^) f + (V„0)' (D^A - KVn) f . (2.48) 
The covariant derivative of a scalar / along the normal n is given by 

V„/ = ^ (/ - N'dj) . (2.49) 

Now we have everything that is needed for the ADM decomposition of the actions of 
CRG. First, in order to shorten the expressions, we denote the decomposition f l2.36p by 

^% = d^(Pd-(P^'^G,,. (2.50) 

Then the action (12.13P is written in terms of the ADM variables: 

Ss = j dtd'x^N \^^^2^1^^-±^ - [D'(f)D^(l) [{D,D, - g^jO'D^ 

-g,j{D'' In N)Dk) % - {K,, - g,,K) Vj% + g^jVnVj%] 
-2 (V„0) D> (AV„ + K,, A) (% + (V„0)' {D'D, - KVn) %] } • (2.51) 

The action (12.140 can be written as 

-g,,{DHnN)Dk) % - {K,, - g,,K) Vj% + g^,VnVJ%] 
-2 ( V„0) A0 (A V„ + KijD^) (% + (V„0)' (D^ A - KWn) ^%T} ■ (2.52) 



We assume that the boundary terms originating from total derivatives vanish when ap- 
propriate boundary conditions are imposed, even though the role of surface integrals is 
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known to be very important in general relativity. Introducing the decomposition fl2.36p 
and the covariant derivative along n fl2.49p explicitly into f l2.5ip would clearly result into 
a very complicated expression. The action f l2.5ip contains time derivatives of the metric 
Qij up to fourth order, time derivatives of the lapse and the shift vector Ni (or of its 
spatial derivative DiNj) up to third order, and time derivatives of (p up to third order. 
Some higher time derivatives of the ADM variables could be removed via integration by 
parts, but with the price of introducing even higher time derivatives of the scalar field (p. 
For the Lagrangian L3 defined by the action f l2.5ip we find 

= 0, (2.53) 



kl 

where g^J^ = ^g§r-- This is clear since the highest time derivative is contained in the term 



a^^^g^{D'(j)Di(f)yVnVnVnK , (2.54) 



where 



V„VnVn-f^ = TTTTzdii + (tcrms with lower order time derivatives) . (2.55) 

Thus q^f cannot be solved in terms of the canonical variables, nor q^^^ defined in terms 
of the initial data from the equations of motion. 

For the Lagrangian L4 defined by the action (12.521) the highest time derivative is 
contained in the term 

- a(DVA0)'(V„V„V„i^^)' (2.56) 

and 

-V9N'-^a{D^<PD^<Pr . (2.57) 



Therefore we expect that g^j^ could be solved in terms of the canonical variables. 

However, since neither of the actions (12.5 ip or (I2.52p depend on iV^^^ , iV-''^^ or (p'^^^ , the 
highest order Hessian matrices 

and , m ' (2-58) 



where qi = (0, A^, A^', gij), / = 1, 2, . . . , 11, have the ranks zero and six everywhere, respec- 
tively. Thus there are eleven and five Lagrangian constraints associated with L3 and L4, 
respectively. As expected, both Lagrangians are singular. The reason for the singularity 
is the presence of gauge symmetry, and its associated (first-class) constraints. 

The actions proposed in Refs. 0, [ij for higher z {z > A) can also be written in the 
ADM form by using the decompositions we have obtained above, (I2.36P and (12.481) in 
particular. Then we see that for the Lagrangian of even z = 2n + 2 and of odd z = 2r;, + 3 
the highest order Hessian matrix is singular, similarly to the two cases discussed above. 
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2.3 Partial gauge fixing 

The constraint (12. 5 p on ensures that the vector d^cf) is timehke everywhere, when 
?7(0) > is assumed. Therefore there exists a preferred fohation of space-time into 
spatial hypersurfaces whose unit normal is given by 

= ^ = ^ [2 59) 

Hence we can write 

a^</.= -v/2f/(0K, d^'(t>=-^2U{<P)n^. (2.60) 

where the unit normal is given in terms of the ADM variables in Eq. (I2.17p . From Eq. 
(I2.60p we see that in this foliation is constant on S^, = 0(t). Thus the constraint (12. 5p 
on reduces to 

-^ + W) = 0. (2.61) 

This implies that the lapse N must be constant on too, = N[t). In order to 
preserve this condition we restrict the symmetry under diffeomorphisms of space-time to 
the symmetry under foliation-preserving diffeomorphisms, given in the infinitesimal form 

as 

5t = f{t), 5x = ^{t,x). (2.62) 

This is the main symmetry group of the HL gravity [7| . In the language of Hof ava's theory 
we would say that both and are projectable — like the lapse A^ is in projectable HL 
gravity. Here we consider (I2.62p as a partial gauge fixing of the diffeomorphism symmetry. 
Now the Eqs. fl2:36|) and (12:481) can be written 

= Uo [K^ - KijK'^ + R) , (2.63) 
(V^V. - (%.V^Vp) / = 2Uo {D'D, - KVn) f , (2.64) 

where (12.71) has again been assumed. Note that is no longer an independent variable, 
since the constraint (I2.6ip implies 

(p{t) = 0(to) + V^o [ dt'Nit') . (2.65) 

Moreover, is not even present in the actions anymore, as we will see next. The actions 
fl23T|) and (l232|) reduce to 

S^ = j dtd'x^N i^^SMEl^^l±^ _ a2Ul {K^ - K.^K'^ + R) 

X {D'Di - KVn) {K^ - Ki,K'^ + R)] . (2-66) 

and 



KijK'^ -K^ + R 



2h? 
—a 



[2Ul {D'D, - KVn) {K^ - Ki,K'^ + R)Y} ■ (2.67) 
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Although these actions have been simphfied considerably, they still contain second order 
time derivatives of the metric gij, and first order time derivatives of the lapse N and the 
shift vector A^*. Indeed, we obtain 



N 



2 ( Kij 



2{K - N'diKj K 

-2 {g^^ - N'dkg'') KuK/ +(R- N^m 



(2.68) 



where the time derivatives of the extrinsic curvature f l2.26p and the inverse metric g'^^ are 
given by 



K.J 



(2.69) 

K = g'^K,, + g'^K,, , g'^ = -g^'g^^m • (2.70) 

We can see that our earlier observations on the singularity of L3 and L4 have not 
changed. For the Lagrangian L3 we obtain 



0. 



For the Lagrangian L4 we obtain 



SgijSgki 



Thus for the Hessian matrices we find 



rank 



6qi6qj 



0, 



rank 



6qi6qj 



6 



(2.71) 



(2.72) 



(2.73) 



and hence we expect ten and four Lagrangian constraints for L3 and L4, respectively. 

Second order time derivatives in the Lagrangian of a theory of gravity are not neces- 
sarily a fatal problem. Indeed we know that actions of the type f{^'^^R) do not suffer from 
the Ostrogradskian instability (for a proper explanation, see the discussion in Ref. [isj]). 
though the scalar curvature (12.251) contains the second order time derivative term 



2V^ (n^K) 



+ (terms with first order time derivatives) 



(2.74) 



In the fi'-^^R) action one can get rid of the second order time derivative by introducing 
additional scalar fields and integrating by parts, although with the price of an extra 
dynamical scalar degree of freedom — the scalaron. In the present case a somewhat 
similar approach can reduce the number of time derivatives in the action. This will be 
discussed next. 



2.4 Further reducing the number of time derivatives 



As discussed in the Appendix of Ref. [5|, in order to reduce the number of derivatives 
in the action, one can introduce four scalar fields (i, ^1, (2 and and write the action 
fl2:T3|l in the form fl: 



d^x^/^g 



2k^ 



a 



[C1C2 + 6 (Ci - (V;.v. - (^.v'Vp) C2) 

+ 6(C2-5^05^0(%.)]} . (2.75) 
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For the action fl2.14p we introduce three scalar field rj, (, ^ in order to write: 

+ e(C-9^0c?>(%^.)]} . (2.76) 

The bad news is that the actions fl2.75p and (I2.76P still contain second order time deriva- 
tives of the scalar fields (2 and (, respectively. 

However, if we employ the preferred foliation discussed in Sec. 12. 3[ we obtain actions 
with only first order derivatives with respect to time. Moreover, the kinetic terms of the 
actions are quadratic in the first time derivatives of the fields. Indeed, we obtain the 
following ADM representations of the actions fl2.75p and fl2.76p : 

dtd^x^N I ^ '^J' ^ - a [C1C2 + 6 (Ci + 2f/o {KVn - D'D,) C2) 

+ 6 (C2 + Uo {g''^'K,jKki -R))]] (2.77) 



and 



S, = j dtd'x^N l^!!!:^!^^ + a[v' + AUoV [KVn - D^D,) ( 

+ UC + Uo {g^'^'K^.Kki -R))]] . (2.78) 
Here we have introduced the so-called De Witt metric 



r^'' = ^ V + - ^? V , (2.79) 



which is the special case a = 1 of the "generalized De Witt metric" 



= I {9'" 9'' + 9'' 9''') - ^9'' 9'' . (2.80) 



When a 7^ 1/3, fl2.80p has the inverse 

'3{a)ijki = - {9ik9ji + 9ii9jk) - 3^ _ ]^ 9ij9kl , g{a)ijkiG{aT'^ = ^^T^'j) ' (2.81) 

Note that the actions fl2.77p and fl2.78p no longer contain time derivatives of the lapse 
and the shift vector N^. This means that the lapse and the shift vector are nondynamical, 
and hence do not propagate, which is also an attractive feature. Obviously the Lagrangians 
of fl2.77p and fl2.78p are still singular due to the remaining gauge symmetry fl2.62p . The 
Hessian matrix for the Lagrangian L3 of fl2.77p is given by 

v^fi-^f^o^l^e-'". 7?^ = -V5"fo4r. (2-82) 



and the rest of the components are zero. The nonvanishing components of the Hessian 
matrix for the Lagrangian L4 of fl2.78p are 



A- A- v^h7^ + "^o^k^^ ' 77I^ = v^"2f/oV. (2.83) 
dgijdgki V2/« / 2iV K^gij ^ 

Thus we expect to have seven and six Lagrangian constraints for L3 and L4, respectively. 
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3 Hamiltonian formalism 



A Hamiltonian formalism for higher derivative theories was first developed by Ostrogradski 
15| . Regular Lagrangians with higher order time derivatives are known to possess negative 



energy degrees of freedom that destabilize solutions of the theory [15|, ll6|. But gauge 
theories are never regular, so Dirac's formalism IT-lij, which is suited for constrained 
(singular) systems, has to be used instead. Generalization of Dirac's formalism to higher 
derivative theories, and a generalization of Ostrogradski's formalism to singular systems, 
were finally achieved in Ref. 20 . 

Lagrangians with higher order time derivatives have also been studied in the context of 
nonlocal theories. In particular see j^, though only regular higher order Lagrangians are 
considered there, whereas all the CRG Lagrangians are singular. Theories where higher 
derivative terms are regarded as corrections to a corresponding lower order theory, and 
especially theories that are truncated perturbative expansions of a nonlocal theory, should 
employ so-called perturbative constraints that ensure that the higher order corrections do 
not change the theory qualitatively (e.g. introduce extra degrees of freedom, lack of lower 
energy bound or violation of unitarity) [igI, 22 . 

GR with a minimal coupling to a scalar (p that is constrained by the same Lagrange 
multiplier constraint (12. 5 p as in CRG has been proposed in Ref. [23 
formalism has been studied in Ref. [2 



Its Hamiltonian 



The Hamiltonian analysis of the actions (I2.5ip and (I2.52p . which contain time deriva- 



tives up to fourth order, could be done using the formalism of Ref. |20|. Likewise, we 
could analyze the simpler actions (I2.66P and (I2.67P or (I2.77P and (12. 78 p . which contain 
second order time derivatives. Here we, however, concentrate on the analysis of the first 
order actions (I2.77P and (12. 78 p . The motivation is to gain some understanding from these 
simphfied CRG actions, because the general ones have proven to be quite difficult to 
analyze. 



3.1 The action 6*3 

First let us define the canonical momenta. Since the action fl2.77l) is independent of the 
time derivatives of iV, iV,, Ci, ^1 and ^2, their canonically conjugated momenta, pn, p\ 
, and , respectively, are the primary constraints: 

Pn^O, p'{x)^0, p^^{x)^0, p^^{x)^0, p^^{x)^0. (3.1) 

The momenta canonically conjugate to gij and (2 are defined by 

V9 {j^^'^'^'^ki - aUo {iig''VnC2 + 66^'^''^^/)) , (3.2) 
-^2aUoiiK. (3.3) 

The Poisson bracket (PB) between the fields and the momenta are postulated in the 
conventional form. The nonvanishing PBs are defined (equal time t is understood): 

{gi,{x),p'^\y)} = 5l5l5\x - y) , {N,p^} = 1 , 

{N,ix),p^iy)} = SiS'ix - y) , {Ci{x),p^,{y)} = " v) , (3-4) 



p-' 

P(2 



SS3 
SC2 
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and similarly for the fields ^i, (2 and ^2 as for (i. Recall that the lapse variable was 
constrained to be projectable, = N{t), due to the constraints (12.611) and = (p(t). 

In order to obtain the Hamiltonian, we perform the Legendre transform in the usual 
manner. First we solve (2 and Qij in terms of the canonical variables (the fields and their 
conjugated momenta): 

/ l-2K^aUo^2 



C2 = N^d,C2 



9ij 



2k^ 



P 



2K^aUo^i 



^l-2K^aUo^2 



Hjkl 



-.kl 



-PC2 



-2K^aUo^2^ 
2K^aUo^i 



-PC, 



where we denote p = Qijp^-' ■ Then we define the total Hamiltonian by 



d X gijp'^ + C2PC2 + ^'nPn + Kp' + AciPCi + hiP^i + -^6^6 " L 



(3.5) 
(3.6) 
(3.7) 

(3.8) 



where Xn, K, A^^, A^^ and A^j are Lagrange multiplier fields. We obtain the Hamiltonian 
by substituting (^-(^ into 



H3 = J d^x {Nn^^ + Ni-Hi + XnPn + Xip' + X^,p^, + X^,p^, + X^,p^,) 



(3.9) 



where we have defined 



^l-2K^aUo^2 
r l + 2K^aUo^2 

+ ^1 2;r^ — 

'Hl = -2D,p^^ +DX2PC,- 



1 /I - 2K'^aUa^2\ 1 /^l - SK^af/oG^ ^ 



2K^aUo^i 



PPC2 + 



6 V 2K'^aUo^i 



R + d [C1C2 + ^1 (Ci - 2f/o/^' AC2) + 6C2] I , 



Pi 



(3.10) 
(3.11) 

Here we denote pij = QikdjiP^^ ■ "^3 are the momentum constraints — three constraints at 
every point on S^. Because of the condition = N{t) the Hamiltonian constraint is the 
integral of l-L^ — a single global constraint. 

The primary constraints (13. ip must be preserved under time evolution generated by 
the total Hamiltonian of the system: 



Pn 


= {pn, Hs} = 


-I 


p' 


= {p\H,} = - 


nji 


P<i 


= {PC,,Hs} = 


-N{ 


PHi 


= {p^-.,H,} = 


-N{ 



-N 



1 - 2K^aUo^2 



P^2 



{p^„H,} = -N{j d'xnl,p^,} 



2K^aUo^i 



pI + Vga (Ci - 2U^D'D,C2) 



(3.12) 
(3.13) 
(3.14) 

(3.15) 



{^9 



2K^aUoi2f 



PijP ■ 



3^ 



(C2 - UoR) 



(3.16) 
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Therefore we impose the following secondary constraints: 

$0 = 1 d'xn', ^ , mx) =ni^O, ^tix) =(2 + 6^0, (3.17) 

<!>Ux) = i (pp,, - ^-^^-^pI) + 3a^f/oe? (Ci - 2f/o/^^A:C2) - , (3.18) 
^6,, ^( ^J I2 1 /^ l-2/t^a^o6 ^' 2^ , (1 - 2/t^a[/o6)% ^ . 

^^(-) = -\P:P'--f J i^^^^^ (C2-f/oi?)-0. 

(3.19) 

Note that in the right-hand side of expressions like (I3.18P we often omit the argument 
X when there is no risk of confusion. Evidently we are dealing with a more complicated 
structure of constraints than in GR or in (projectable) HL gravity. We will see that the 
constraints $3 (a;), n = 4, 5, 6 are second-class and they enable one to define the auxiliary 
fields ^i,Ci5^2 in terms of the dynamical fields gijX2 and their canonically conjugated 
momenta p^^ , p^^ . 

It is convenient to introduce a global smeared version of the momentum constraints 

^l{x^)= f d'xxiHl (3.20) 



where (^ = 1, 2, 3) are arbitrary functions on which vanish rapidly enough at infin- 
ity. As the name suggests, the momentum constraints generate the infinitesimal spatial 
diffeomorphisms of the dynamical variables: 

{MiXk),gij} = -x^dkgij - gikdjX^ - gjkdiX^ , 
{$f(x.),C2} = -x*5,C2, 

{^l{Xi),P^,] = -diX'PC2 - X%PC2 ■ (3.21) 

Its PBs with the nondynamical fields and their canonically conjugated momenta vanish. 
However, when ^i,Ci,^2 are solved in terms of the dynamical variables, they behave as 
scalar fields, i.e. like (2 in Eq. fl3.2ip . The momenta are tensor or scalar densities of 
weight —1 under spatial diffeomorphism, while the fields are regular tensors or scalars. 

In order to check that the secondary constraints $3 (/ = 0, 1, 2, . . . , 6) are preserved 
under time evolution, we need to evaluate their PBs with every constraint. First consider 
the momentum constraint (13.201) . Its PB with the primary constraints (13.11) vanish. Since 
^{3 is a scalar density of weight —1 under spatial diffeomorphism, its integral vanished 

{^|ix^),^l} = 0■ (3.22) 
The PB of the momentum constraint with itself forms the Lie algebra 

{^liXi), $f (^,)} = ^lix'd,^^ - i^'^,x^) ^ , (3.23) 

Then consider the PBs with the rest of the secondary constraints in Eqs (I3.17p - (l3.19p . We 
see that the constraints ^^{x),n = 4,5,6 have been defined to be scalars under spatial 
diffeomorphism. 



*For a scalar density of weight w we find: {^f (Xi)i/ d'^xifj^^} = —{w + 1)/ d'^xx^diipti 
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The Hamiltonian constraint $3 has vanishing PBs with and p*. The PBs with the 
rest of the primary constraints were calculated in Eqs. f l3.14p -( l3.16p . and they vanish due 
to the secondary constraints fl3.17p - fl3.19l) . Since the Hamiltonian constraint is global, its 
PB commutes with itself 

{*3,*3}=0. (3.24) 
Next consider the PBs with the secondary constraints $3 (as), n = 4, 5, 6: 



P 



1 - 2fi;2af/o5 



(3.25) 



and the PBs of the last two secondary constraints turn out to be rather complicated 
expressions: 



3k^ 



g3/2 I _ 2K^aUoi2 



PijP 
2 



1 2 1 / 1 - 2K2af/o^2 



nP + n 



Pi 



PC2 + 



1 r 1 + 2K2a[/o6 



4k2 



PPC2 



R 



a 



+ ii ^Ci - 2UoD'D,C2 + ;t6C2 + Uo (^,^1^X2 - 2D^A6) 



PC2 



2k^ 

1 - 2K2af/o6 



2p„- -^{2p 



1 - 2K^aUo^2_ 
2K^aUoii 



-PCs 



1 



'3af/o6 



_ 1 - 2K2af/o^2 



(3.26) 



and 



nP -n 



1 fl- 2K'^aUQ^2 



6 V 2K^aUo^i 



PC2 



+ 



1 r l + 2K^aUo^2 (2 



V9 



2^2 



-Rp - Rijp'^ ) - a^2C2P 



2 ( P'' - -P^?*^' ) (l^./^,6 - 2D^,^,D^)C2) - 26 ( 



+ -{Ci + ^2-2UoD'DS) 
{l-2^^aUo^2y 



P 



1 - 2K'^aUo^2 
2K^aUo^i 

1 - 2K'^aUo^2 
2K^aUQ^i 



PC2 



-PCs 



2k' 



2p^_ 



9ij 



l-2ifaUg2 
2.2a% 



(3.27) 
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The nonvanishing PBs between ^^{x),n = 4,5,6 and the primary constraints (13. ip 

are 

{<l>t{x),pM} = S'i^ - y) , mx),pUy)} = Sa'UoeiS'ix - y) , 
{^\{x),v^,{y)} = -g^/c,S%^ - y) , (3-28) 

Now we can consider the stabihty of the secondary constraints. First we obtain that 
the Hamiltonian constraint and the momentum constraints are preserved under time evo- 
lution: 

= {^l H,} = ! d^x^g U^a^l + A,, -1-,$^ + A, - , 

J V otttyo'Si (1 - 2K^aUQ^2> J 

^{X^) = {<^|{X^),H,} = {$f(x.),$f(iV.)} = ^lix'd.N, - NW,x^) ^0. (3.29) 
The rest of the secondary constraints evolve under time as follows: 

mx) = {^l{x), H,} = N{^l{x), $°} + N^^l + Xc,3a'Uoei 

mx) = mx), H,} = N^x), $°} + N^d.^l + A,,i^ {^^^U^JpI 

where Eqs. (I3.26P and (I3.27P are understood. We see that imposing ^^{x) ^ 0, n = 4, 5, 6 
in Eq. (I3.30p fixes the Lagrange multipliers A^^, A^^ and X^^- Thus no more constraints 
are required to fulfill the consistency conditions for the secondary constraints. Hence we 
have obtained a stable constraint surface in the phase space. 

As noted above, we may use the constraints $3 (a;), n = 4, 5, 6 to eliminate the auxiliary 
fields Ci, ^1 and ^2- To this end, we can set the second-class constraints p^^, p^^, p^^ 
$3(a?),n = 4,5,6 to vanish strongly by introducing the Dirac bracket, and then solve 
the auxiliary fields in terms of the dynamical variables from Eqs. fl3.17p - fl3.19p . This 
eliminates the fields Ci,C,i,C,2 and their canonically conjugated momenta from the phase 
space of the system. 

Let us then count the number of physical degrees of freedom (physical d.o.f.) by using 
Dirac 's formula: 

^(physical d.o.f.) = - [^^(canonical variables) — 2 x ^^(first-class constraints) 

(second-class constraints)] . (3.31) 



16 



For nonpropagating zero modes we have 28 canonical variables (A^, Ni, gij, ^i, ^i, ^2, C2, and 
their conjugated momenta), 8 first-class constraints {p]\f,p\^'^,'Hl), and 6 second-class 
constraints (^(^-^,^5^, j9^2) "^3) "^35 "^s)- This yields 

# (nonpropagating physical d.o.f.) = ^^'^^ — 16 — 6) = 3 (3.32) 

for the (a;-independent) zero modes. For propagating modes we have 26 canonical variables 
{Ni, Qij, ^1, (i, ^2, C25 and their conjugated momenta), 6 first-class constraints {p\ T-Ll), and 
6 second-class constraints (p^^ , , , $3, $3). Thus for propagating modes we obtain 

^(propagating physical d.o.f.) = 2 (^^ — 12 — 6) = 4 . (3.33) 

For comparison, there are only two physical degrees of freedom in GR. As another com- 
parison, our analysis shows that CRG with the condition = 0(t) has one more physical 
degree of freedom than projectable HL gravity, which has 2 zero modes and 3 propagating 
modes. Interestingly, the number of physical modes is exactly the same as in the modified 
F{R) HL gravity Q, 3 lll| . One extra physical degree of freedom has its origin in the 



higher order time derivatives present in the CRG action. The other extra propagating 
mode is caused by the projectability condition similarly as in HL gravity. Such extra 
degrees of freedom can be problematic since they may generate extra (long range) forces 
that are not in agreement with observations. One may be able to bring the number of 
physical degrees of freedom closer to that of GR by introducing some extra gauge symme- 
try, which generates some new constraints. Another possible way to deal with the extra 
propagating mode is to make the scalar field (2 massive. 



3.2 The action S4 

Since the action (12.781) is independent of the time derivatives of A^, A^*, rj and ^, their 
canonically conjugated momenta, p^, p\ Prj and p^ respectively, are the primary con- 
straints: 

PN-0, p\x)^0, Pr,ix)^0, ps:ix)^0. (3.34) 
The momenta conjugate to and ( are defined by 



6S. 



PC 



5C 



6S4 



The PB is postulated similarly as in Eq. (13. 4p . 

Solving ( and Qij in terms of the canonical variables gives 



N 



9ij 



1 l + 2K^aUo^ 



1 + 2K^aUo^ 



2NKi, + 2Z}(,Ar,.) 



2k2 



Hjkl 



1,1 g^^ { 1 + 2K''aUQi 



PC 



(3.35) 
(3.36) 



(3.37) 
(3.38) 
(3.39) 
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The total Hamiltonian is obtained similarly as before: 



(3.40) 



where the Hamiltonian and momentum constraints are defined by 



-1/0 

7x4 



2k' 



2^- 

•HI = -2D,p'^ + DXpc 



PijP' 



:P + 



1 / 1 + 2K''aUoi 



3' 3 V An'^aUor] 



PPC + 



1 /l + 2/t2a[/o^ 



6 \ AK'^aUor] 



(3.41) 
(3.42) 



The primary constraints fl3.34p must be preserved under the time evolution generated 
by the total Hamiltonian H4 of the system: 



Pn = {pn, Hi} 



4 1 



P 



{p\Hi} = -ni, 

Pr, = {p^,Hi} = -N{ j d^Xnlpr,} 

' 1 1 ( \^2K^cdJ^i 



N 



p^ = {p^,H,} = -N{ j d'xnlp^} 



pI + ^2a {t] - 2UoD'D,C) 



(3.43) 
(3.44) 

(3.45) 



N 



V9 



PijP ■ 



pI 



2AaUori'' 



(C - UoR) } . (3.46) 



Therefore we impose the following secondary constraints: 



d-^xHl ^ , 



^l{x) 



PijP 



1 A +2/t2a[/o^ 



3^ 6 \^ An'^aUof] 



Pi 



[1 + 2K^aUoO^ 



(3.47) 
(3.48) 

{(-UoR)^0. 

(3.49) 



The constraint structure is quite similar to the case z = 3 we discussed above. The 
only major difference is that there is one auxiliary field less in the case z = 4 than in 
the case z = 3, and therefore there are two second-class constraints less than in the case 
z = 3 — one primary constraint and one secondary constraint. The consistency conditions 
$4(3:;) ^ 0,n = 4,5 can be satisfied by fixing the Lagrange multipliers and Ag. Thus 
we again obtain a stable constraint surface. Once imposed, the second-class constraints 
fl3.48p and fl3.49p define the auxiliary fields 77 and ^ in terms of the dynamical variables. 
The number of physical degrees of freedom is the same as in the case z = 3. 
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4 Conclusion and discussion 



Let us summarize our approach and results. Following the proposal of Ref. we consid- 
ered that the exotic perfect fluid of CRG is generated by a scalar field (f) that is constrained 
by Eq. (12. 5p . We obtained the ADM representations of such CRG actions for the sup- 
posedly power-counting renormalizable 2; = 3 and super-renormalizable z = A cases. The 
corresponding Lagrangians were found to contain time derivatives of the dynamical fields 
up to fourth order. In order to obtain Lagrangians with kinetic parts quadratic in the 
first order time derivatives, we introduced some additional scalar fields and took advan- 
tage of the constraint (12. 5p by choosing d^cj) to be parallel to the normal of the spatial 
hypersurfaces. The Hamiltonian formalism for the space-time decomposed actions was 
developed. We showed that both actions have a consistent set of constraints. All but 
one of the additional scalar fields turn out to be auxiliary, which can be eliminated by 
imposing the second-class constraints of the systems. The number of physical degrees of 
freedom was found to be 3 for zero modes and 4 for propagating modes. Thus we can 
conclude that the advantage of retaining general covariance and local Lorentz invariance 
in the formulation of CRG comes with the price of an extra physical degree of freedom 
— compared to HL gravity. This is a consequence of the higher order nature of the CRG 
action. Compared to GR we found two extra propagating physical degrees of freedom. 
One may be able to reduce the number of physical degrees of freedom closer to that of 
GR by introducing some extra gauge symmetry or by making some of the fields massive. 
These points, however, were not pursued further in this paper. 

It would be interesting to use the Pons formalism ||20|] for the Hamiltonian analysis of 
CGR in order to gain deeper understanding of these higher order derivative theories, and 
with the hope that a more conclusive statement on their nature could be drawn. It would 
enable one to analyze CRG without imposing the projectability condition on the scalar 
field 0, and consequently on the lapse N . Another possible prospect is the application of 
perturbative constraints to CRG. This could enable one to eliminate the possible problems 
caused by the higher order time derivatives. It would also be interesting to see whether 
the "covariant vector gravity" proposed in Ref. ^ could provide even more fruitful results 
compared to the present scalar formulation, since unlike 9^0, a vector field does not 
involve time derivatives. 
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